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Abstract 



We argue that in quark matter at high densities, the color magnetic field 
remains unscreened and leads to the phenomenon of color superconductiv- 
ity. Using the renormalization group near the Fermi surface, we find that 
the long-range nature of the magnetic interaction changes the asymptotic 
behavior of the gap A at large chemical potential \x qualitatively. We find 
A ~ [ig~ 5 e~x.p(— • |), where g is the small gauge coupling. We discuss 
the possibility of breaking rotational symmetry by the formation of a conden- 
sate with nonzero angular momentum, as well as interesting parallels to some 
condensed matter systems with long-range forces. 
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I. INTRODUCTION 



The study of QCD at finite density has a long history. The suggestion that at high 
densities hadronic matter becomes quark matter IJ was made almost immediately after 
the discovery of asymptotic freedom 0. It has been known for almost as long that at 
very high baryon densities, where perturbative QCD can be applied, quark matter is in 
a BCS-type superconducting state ||. More recently, various groups 0,|5|| have modeled 
quark matter at intermediate densities using phenomenological four-fermion interactions and 
found that condensates form with sizable gaps of order 100 MeV. These studies reveal a rich 
phenomenology ||J-§], including an interesting phase diagram |6[], new symmetry breaking 
schemes like color-flavor locking Jf], which might be also relevant for the study of nuclear 
matter in the nuclear superfluid phase, if the latter is continuously connected to the quark 
matter phase ||. 

Naturally, a systematic calculation in perturbation theory is possible only in the high- 
density regime, where the chemical potential fi is much larger than Aqcd and the strong 
coupling at the scale of the Fermi energy is small. However, even in this case reliable 
calculations have so far been hampered by the inability to take into account the long-range 
nature of the color magnetic force, which dominates the interaction between the quarks at 
large distances. Instead of tackling the problem of condensate formation by a long-range 
interaction, many treatments rely on the assumption that a magnetic mass of order g 2 \i 
is somehow generated. However, in reality, the magnetic field remains unscreened in the 
absence of superconductivity itself. Ordinary BCS theory cannot be directly applied due 
to the IR divergence from the exchange of the soft magnetic gluons. Therefore, even the 
asymptotic, weak-coupling behavior of the gap has not been found. 

In this paper we present what we believe to be the correct estimate for the value of the 
BCS gap at asymptotically high densities. Our approach is based on the renormalization 
group around the Fermi surface [§,|l(J> properly modified to take into account the long-range 
magnetic interaction. We find that the gap is proportional to exp(— c/g), c = ttj-Q This 
behavior is different from the naive expectation from BCS theory, which predicts exp(— c/g 2 ). 
The fact that the suppression is parametrically much milder means a larger value of the gap 
at high densities, and potentially could also indicate that at intermediate densities the gap 
and the critical temperature may be larger than previously estimated. The latter may 
substantially enhance the chance of forming a color superconductor in heavy-ion collisions. 
We also show that the asymptotic behavior of the gap does not depend on the angular 
momentum of the Cooper pairs, and comment on the possibility of breaking rotational 
symmetry by the formation of a condensate with a nonzero angular momentum. 

We will first review the renormalization group approach to the BCS theory (Sec. [TT]), 
then describe the trouble caused by the long-range magnetic interaction (Sec. |IIJ). In Sec. 
P"V| we describe our resolution of this problem. We then make final remarks in Sec. |V|. The 
Appendices contain various technical details, including a treatment of Eliashberg theory. 



1 This agrees with a comment in a recent paper by Pisarski and Rischke [11] that A ~ exp(— c/g) 
for some constant c. 
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II. THE RENORMALIZATION GROUP NEAR THE FERMI SURFACE 



An elegant method to see the formation of the BCS superconducting state is the renor- 
malization group (RG) near the Fermi surface 0. This approach has been applied to the 
case of quark matter by Evans, Hsu and Schwetz, and Schafer and Wilczek ||10|| . These 
treatments apply when there exists a nonzero magnetic mass, screening the color magnetic 
interaction. For completeness, we give here an elementary overview of the approach, in the 
spirit of non-relativistic quantum mechanics. Readers who need a more formal, rigorous 
treatment should consult Ref. [10]. In this section we will follow Refs. and consider 



a theory of quarks interacting via a local four-fermion interaction. The aim is to give an 
overview of the conventional RG approach before tackling the more difficult problem of a 
long-range interaction. Keeping in mind that one-gluon exchange conserves helicity, we will 
for simplicity consider only left-handed fermions. 

Let us imagine a Fermi gas of massless quarks with a chemical potential fi. In the ideal 
gas approximation, all energy levels below the Fermi surface |p| = fi are filled. The energy 
will be measured relative to the Fermi surface, so we introduce e p = p — fi. 

The RG procedure works as follows. At any given step, the effective theory contains 
only the fermion degrees of freedom located in a thin shell surrounding the Fermi surface. 
These fermions have |e p | < 5. All others have been integrated out. The only relevant 
interaction between the fermions is the scattering of pairs with opposite momenta p|,|T0 



Let us introduce the scattering amplitude from a pair with momenta (p, — p) to another 
pair with momenta (k, — k), 

f(9) = /(p, k) = T(p, -p - k, -k) 

To avoid complications with statistics, we will assume that the two particles are of different 
flavors. Near the Fermi surface, the amplitude depends only on the angle 9 between p and 
k. A positive / corresponds to a repulsive interaction, and a negative / means attraction 
between particles with opposite momenta. 

In the spirit of the Wilson renormalization group, let us now integrate out all fermion 
states with e _1 <5 < |e p | < 5. According to quantum mechanics, the scattering through virtual 
states in this region gives a correction to the scattering amplitude. To account for these 
virtual processes, we need to correct the scattering amplitude. Thus / obtains a correction, 

/(p,k)^/(p,k)-E r(p '' p ;' )7 : ( ^ k -' k) a) 

where the sum is over all intermediate states i belonging to the sector of the theory that 
has been integrated out. The virtual state i may have a different energy than the initial 
energy 2e p . We assume that the initial and final particles are almost exactly located at the 
Fermi surface, so e p = e k = 0. What could be the states il To answer this question one 
notices that the scattering through an intermediate state can be of two types: 

1. The pair (p, — p) can scatter to an intermediate pair (p', — p'), which then goes to 
(k, — k). In this case, the intermediate state % is that with two particle excitations with 
momenta ±p'. The Pauli principle requires that p' is located above the Fermi surface. This 
state has = 2e p / and T(p, — p —>•«) = /(p, p'), T(i — > k, — k) = /(p', k). 
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2. Alternatively, first a pair of particles inside the Fermi sea with momenta (p', — p') can 
scatter to make the final pair (k, — k), and then the initial pair (p, — p) scatters to fill the 
holes vacated by the pair (p', — p') in the Fermi sphere. In this case, the intermediate state 
% consists of six elementary excitations: four particles with momenta ±p and ±k, and two 
holes with momenta ±p' located below the Fermi surface, p' < \i. In this case, Ei = — 2e p /, 

r( Pl - P -> i) = f( P ', k), T(i -> k, -k) = /( P , P '). 

It is clear now that Eq. (|I|) becomes 

/(p,k)^/(p,k)-/ / < p -p; )/(p, - k) ( 2) 

Jp' ^Fp'l 

where the integration is over all p' satisfying e~ 1 5 < \p — fi\ < 5. The integral over |p'| can 
be taken, and Eq. (0) reads, 



/(p, k) -> /(p, k) - £- 2 j ^/(p, p')/(p', k) 



where the integration is over all directions of p'. Repeating the RG procedure many times, 
one finds that / evolves according to the RG equation, 

where t = — In S goes to +oo as one approaches the Fermi surface. Eq. ([|) describes the RG 
evolution of the scattering amplitude on the Fermi surface. 

It is convenient to expand the scattering amplitude over partial waves, 

oo 

=£(2/ + l)/ifl(cos0) 

1=0 

or, inversely, fi = \ j^dO sin 6 P;(cos 6)f{6). Using a well known property of the Legendre 
polynomials, we find that the partial-wave amplitudes // evolve independently, 



dfi fi- 



2 



dt 2tt 2 



ft (4) 



The solution to Eq. (01) is 



fl(t) 



i+&yi(o)t 



We see that if at t — all // > 0, which means that the interaction is repulsive in all 
channels, then the four-fermion interaction vanishes at the Fermi surface. However, if one of 
fi(0) is negative, it will develop a singularity (the Landau pole) at t — — ■ The Landau 
pole is reached first by the channel having the largest negative fi(0). This singularity is 
nothing but the manifestation of the BCS instability of the Fermi surface with respect to 
any attractive interaction. The BCS gap is proportional to the energy scale at which the 
Landau pole is reached, 
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2tt 2 



^ 2 /(0) 



A ~ exp 

Let us reproduce some results obtained in Ref. fL0| . Let us take the interaction in the 
form G°{ip r )°ip) 2 + G 1 ^^ 1 ^) 2 , where G° and G l are two independent constants. The tree- 
level scattering amplitude p, — p — > k, — k between two left-handed particles arising from 
this interaction is 



f{9) = 2 



G°cos 2 --G i fcos 2 - + 2sin 2 - 



2 V 2 2 



(5) 



which contains only the s-wave and p-wave terms, with fo = G° — 3G l and f\ = + G 1 
Therefore, G° and G 1 can be said to run according to the following equations, 



^(G -3G i ) = -^(G t, - iG f 



1(G° + G l ) = -^(G° + G l 
at 6tt 2 



which constitute a subset of the equations found in Ref. [10 . 

As a toy model mimicking the real one-gluon exchange, one can follow Ref. [0 and 
take the interaction of the form — -^(ip'y^ip) 2 . Here A should be thought of as the typical 

2 _ 

momentum of the exchanged gluon, and \ is the effective coupling in the color 3 channel 
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T0|1 j where the superconductivity effect is usually strongest. This corresponds to G° 



2 

-G l = —^2- The interaction is most attractive in the s-channel and a BCS state is formed 



with the gap 



, 3ttA 2 , 



The gap has the e c ^° 2 dependence on the coupling g. This parametric dependence of the 
gap has been obtained in Ref. |10j and is the same as that obtained in variational or mean 



field treatments of models with a four-fermion interaction between quarks We will see 

that in the true theory, the dependence at asymptotically high densities (g — > 0) is different 
from that obtained in these toy models. 

III. THE PROBLEM OF THE UNSCREENED MAGNETIC FIELD 

Let us try to naively apply the method developed in Sec. || to high-density QCD. At the 
lowest order, the quarks interact by exchanging one gluon. The gluon propagator, which is 
1/q 2 in vacuum, is modified by screening effects. The static color electric field is subjected to 
Debye screening at the distance scale m^ 1 , where r«D ~ g/J,, which can be seen by resumming 
bubble diagrams in the gluon propagator. In the magnetic sector, the same resummation 
yields a magnetic gluon propagator, 

D{qo ' 9) = tfi + 'm* M (?) 
q -+- 2 m D~ 
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in the regime go g The term | comes from the Landau damping. In the static 

limit g = 0, the magnetic field is not screened. If go 7^ 0, the field is said to be "dynamically 
screened" on the scale g ~ rn^ 3 q^ 3 . 

Before going further, let us make an important comment on a confusion in the literature. 
This confusion originates from the similarity between the high-density and high-temperature 
gauge theories. The similarity can be seen in the Debye screening, which occurs at the scale 
gT at high temperatures and gfi at high densities. At high temperatures, the magnetic field 
is not screened at the one-loop level, but develops a non-perturbative screening of order 
g 2 T. This sometimes leads to an unjustified assumption that the magnetic field develops a 
magnetic mass of order g 2 [i in high-density QCD. 

To see why the analogy between high temperatures and high densities breaks down on 
the question of the magnetic screening, let us review three standard ways to interpret the 
emergence of the magnetic mass in hot gauge theories. The first argument is dimensional 
reduction: the static large-distance behavior of a gauge theory at high temperatures is the 
same as of an effective Euclidean three-dimensional gauge theory. The latter is confined 
on the scale g 2 T, which means that the magnetic field in the original theory should also 
be confined at this scale. The second argument is that the high-temperature perturbation 
theory is IR divergent for momenta < g 2 T, and only a magnetic mass of this order could 
make the perturbation theory finite again. The third way is to notice that, due to the 
Bose enhancement, the thermal fluctuations of the gauge field become so large at the scale 
g 2 T that they are fully non-linear. The last argument does not necessarily imply magnetic 
screening; it just shows that the existence of the latter at the scale g 2 T does not contradict 
perturbative results, since at this scale the physics is non-perturbative. 

None of these three arguments can be carried over the case of high densities. First, there 
is no dimensional reduction for gauge theories with a finite chemical potential. Second, the 
perturbation theory for the long-range magnetic field is infrared finite. Indeed, in vacuum, 
the infrared divergences of one-loop graphs come from integrals like / d A qD 2 (q), where D(q) ~ 
q~ 2 is the gluon propagator. At finite fi, the gluon propagator is modified as in Eq. (|7|). Now 
go effectively scales like g 3 instead of g, and by a simple power counting one sees that the 
integral is finite in the IR. Therefore perturbation theory does not break down for momenta 
of order g 2 /j,, or g n fi with any n, and there is no reason to expect non-perturbative effects 
proportional to any finite power of g^\ The third argument explicitly relies on the large Bose 
enhancement that takes place only at finite temperatures; this argument clearly does not 
work at finite densities. Therefore, the magnetic interaction is not screened at the scale g 2 \i 
as is occasionally assumed. 

Now let us return to the our problem and try to apply the formalism developed in Sec. 
[TI] to the interaction mediated by one-gluon exchange. We will see immediately that we 
have serious trouble with the very soft gluons. Indeed, on the Fermi surface, the tree-level 
small-angle {6 <C 1) scattering amplitude, due to one-gluon exchange, is 



2 The author thanks S.Yu. Khlebnikov for pointing out this argument. Similar conclusion has been 



reached by Pisarski and Rischke [11]. 
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ftree{0) ~ f~ ( 2Q2 jf ^ 2^2) (^) 

3 \fi z 9 z + mf ) fi Z 9 2 / 

The two contributions in the RHS come from the electric and the magnetic interaction, 
respectively (again, the factor | comes from considering only the 3 channel.) All partial 
amplitudes diverge logarithmically. For example, 

f = - [ d9sm9P l (cos9)f(9)*-?-ln-^- (9) 

IJmm/M 

where g m i n is the smallest allowed momentum exchange that one has to put in by hand to 
make /o finite. Clearly, this IR problem renders the conventional RG formalism unusable. 

Let us also warn against what might seem to be a natural solution to this IR problem. If 
one assumes that a quark condensate is eventually formed, the magnetic field is screened by 
the Meissner effect. One may be tempted to take the inverse London penetration length g A 
as the cutoff in Eq. @, and use the computed value of f to find the gap, thus obtaining a 
self-consistency condition for A. In this way one does obtain a gap of order e~ c ^ 9 , where c is 
some constant. However, this approach is flawed, and gives the wrong value of c, because it 
entirely neglects the screening effect of Landau damping which turns out to be much stronger 
than the Meissner effect. To see this, note that the condensate smears out the Fermi surface 
over a scale A, which is the natural energy scale of quasiparticles near the Fermi surface. 
The gluons that such excitations exchange have go ~ A. On these frequencies, the dynamical 
magnetic screening happens already at the scale q ~ m^A 1 / 3 , which is much larger than A. 
Therefore, the source of the IR cutoff should be the dynamical magnetic screening, which 
already takes place in the normal phase, rather than the Meissner effect. 

Let us now turn to the central part of this paper, where we will find the correct RG 
treatment of the theory with the color magnetic interaction. 



IV. RENORMALIZATION GROUP FOR MAGNETIC INTERACTIONS 

Let us repeat the RG procedure described in Sec. ||. At any given step in the RG 
evolution, one keeps only fermion modes having energies smaller than 5. The start of the 
RG evolution, t = 0, will be taken at 5 ~ mr>, and the evolution stops when 5 is of order 
the gap, so typically 5 <C m^. 

At tree level, the fermions interact via one-gluon exchange, characterized by the momen- 
tum of the gluon (q , q). Since all fermions have energy less than 5, the energy of the gluon 
go is naturally of order or less than S, while the momentum exchange q can be anywhere 
between and 2/x. 

Let us divide the four-fermion interaction that arises from the one-gluon exchange into 
"instantaneous" and "non-instantaneous" parts. The instantaneous interaction is mediated 
by the gluons that have momenta q > q$ = m^d 1 ^. The Landau damping for these gluons 
is negligible, m^^- < q 2 . The gluon propagator, which is now simply q~ 2 , does not depend 
on q , which means that the four-fermion interaction they mediate can be considered as 
instantaneous. This part of the interaction is of the familiar type and will be treated in 
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the conventional way. In particular, one can characterize this part by the partial-wave 
amplitudes //. For q < q$, the Landau damping can no longer be neglected. This part of 
the interaction has a considerable temporal retardation and should be treated separately. 

Now let us reduce 5 by a factor of 1/e by integrating out fermion degrees of freedom 
with energy between e~ l 5 and 5. During this process the following will occur: 

1. The partial- wave amplitudes /; obtain the conventional renormalization, as written 
in Eq. (g). 

2. One could ask if the non-instantaneous coupling is renormalized during this integra- 
tion. To answer this question, one should compute the correction to the non- instantaneous 
interaction that comes from integrating out the fermion degrees of freedom. In Appendix [A] 
we explicitly estimate the corresponding one-loop diagrams and show that the non-instan- 
taneous part of the interaction does not get renormalized. 

3. Most importantly, and what makes our RG distinctive, part of the non-instantaneous 
interaction becomes instantaneous. Specifically, the gluon exchange with q lying in the 
interval (e~ 1//3 g<5, q$), which was formerly treated as non-instantaneous, now becomes a part 
of the instantaneous interaction and contributes to f]. Simply speaking, our criterion of 
what to consider as instantaneous has became more inclusive, since we are now looking at 
a smaller energy scale, corresponding to a larger time scale. 

How much of the non-instantaneous part of the interaction transfers to the instantaneous 
part during one step of the RG? According to Eq. (H) and the non-renormalization of the 
non- instantaneous interaction, the increment in f(9) has the form, 

and vanishes outside this window of 9. For definiteness, let us concentrate our attention to 
the s-wave amplitude /q. This amplitude obtains a constant additive contribution from the 
soft sector at each RG step, 

A/ = - / d9 sinO A f (9): 

Therefore, the RG group equation for /q now becomes, 



d_ f 9 2 M 2 
dt Jo V 2tt 2 ' 



TtJ® - -7^2 ~ TT^fo (11) 



The second term in the RHS is the familiar term that gives rise to the BCS effect for short- 
range interactions. What is new is the first term, which takes into account the fact that softer 
and softer gluon exchanges contribute to f). The non-instantaneous part of the interaction 
can be considered as an infinite pool, which continuously replenishes the instantaneous part 
during the RG evolution. Clearly, this should speed up the approach to the Landau pole. 

To secure a solution we also need to specify an initial condition on / . Recall that t = 
corresponds to 5 ~ jtid, from Eq. (|$D one finds, to the leading logarithm, 

/o(0) = -|4 ln - (12) 
3/i 2 g 
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The solution to Eq. ( |TT| ) with the initial condition (|12|) is 



Jolt) = — ^~ tan 



9 ^t + ein 1 



3y/2ir\ g 



The coupling / hits the Landau pole when the argument of the tangent is equal ?. This 
happens when 

3vr 2 , 1 
t = —=- - 6 In - 

The Fermi liquid description, thus, breaks down at the energy scale 

A ~ m D e~* ~ fxg~ 5 exp ( 13 ) 



which will be interpreted as the scale of the gap. Notice that the gap is proportional to 
e~ c / ff , which is parametrically larger than the naive estimate e~ c ^ 9 at small g. The reason 
for this enhancement is obviously the singularity of the magnetic interaction. 

Strictly speaking, the RG calculation does not tell us that A is the value of the gap. 
In fact, the RG merely indicates that the normal Fermi liquid behavior breaks down at 
the scale of A. To confirm that A is the gap, one needs to use some alternative approach. 
In Appendix |B|, by making use of the Eliashberg equation, borrowed from the physics of 
electron-phonon systems, we verify that the gap A is in fact proportional to exp(— -^). 

V. REMARKS AND CONCLUSION 

We have found that the ground state of the system of quarks, interacting via one-gluon 
exchange, is basically a BCS-type superconductor, and found the weak-coupling behavior 
of the gap using the RG approach, appropriate in the presence of a long-range magnetic 
interaction. We found A ~ g~ 5 e~ c ^ 9 , c = not A ~ e~ c ^ g as in conventional BCS theory. 
Thus, at least in the g — > limit, the superconductivity gap is larger than in all previous 
estimates P~|^,pi|. Here let us make a few remarks on our calculation. 

Wave function renormalization and non- Fermi-liquid behaviors. A similar problem of 
fermions interacting via an unscreened U(l) magnetic field has also generated considerable 
interest in condensed matter physics. The fermions may be the valence electrons in metals 



12|1 , in which case the U(l) field is the magnetic component of real electromagnetism, or 



some effective degrees of freedom in low-dimension strongly correlated electron systems |13 
where the U(l) gauge interaction could be generated as an effective interaction. In metals, 
the magnetic interaction is repulsive for a pair located on the opposite sides of the Fermi 
surface, so it does not lead to BCS superconductivity. However, interesting phenomena may 
still arise from this repulsive interaction. In particular, in has been shown [12[] that the weak 



magnetic interaction leads to the breakdown of the Landau theory of Fermi liquid, typically 
at extremely low temperatures. One could ask whether the effects leading to this non-Fermi 
liquid behavior would modify our calculations. 
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The basic idea is that the fermion wave function gets a large renormalization near the 
Fermi surface from the magnetic interaction. To one-loop level, the renormalization of the 



wave function is [12 



Z- l {q ) = 1 + const -g 2 ln— (14) 

If one goes arbitrarily close to the Fermi surface, the wave function renormalization Z tends 
to 0, which means that the discontinuity of the distribution function at the Fermi surface 
disappears, thus signaling a deviation from Landau's Fermi liquid theory. However, in our 
case, the BCS effect is already essential at go ~ A. Recalling that A ~ e~ c ^ 9 , we find from 
Eq. (|TJj) that Z remains close to 1, Z — 1 + 0(g). Therefore, one can safely ignore the 
renormalization of the wave function, and there is no chance for non-Fermi-liquid behavior 
other than BCS-type superconductivity to manifest itself. 

We note here that it was suggested that an attractive magnetic-type interaction may 
arise in various situations in condensed matter physics, for example in double-layer electron 
systems. A BCS gap may emerge in such systems. Due to a stronger singularity of the 
magnetic interaction in 2d, the gap is found to be proportional to a power of the coupling 
constant, rather than being exponential JT4j . 



The very soft magnetic gluons. In our RG approach, even after the final step of the RG 
evolution, the interaction still contains a non- instantaneous sector, which is carried by the 
magnetic gluons with energy go ~ A and momentum q < m^A 1 / 3 . The BCS effect is due 
to the instantaneous interaction, but one may ask whether the remaining non- instantaneous 
interaction could destroy the BCS state. Here we give an (admittedly crude) argument as 
to why this cannot happen. 

The magnetic modes mediating the non-instantaneous interaction have their intrinsic 
time scale g^ 1 ^ A -1 . Therefore, these modes can be considered as static during the typical 
time scale of the system, A" 1 . The question is now whether this random, static magnetic 
field could destroy the superconducting state. Let us estimate the total strength of the 
quantum fluctuations of the magnetic field. It is roughly 

is 1 
B 2 ~ [d q q 2 \A(cO\ 2 ~ [dq dqq 2 - ~ m 2 D A 2 

so the typical value of the fluctuating magnetic field is mj)A ~ gfiA. This should be 
compared with the critical magnetic field, which is of order /xA. We conclude that at weak 
coupling, the almost static quantum fluctuations of the magnetic field are too small to 
destroy the superconducting state. 

The possibility of breaking rotational symmetry. In our treatment of the RG equation, 
we have concentrated our attention on the s-wave coupling fo. However, we could equally 
consider higher partial waves. Let us take arbitrary I and try to rewrite Eq. ([TTj) for //. 
According to Eq. @, the f 2 part in the evolution of // has the same coefficient as for f Q . 
The constant part depends on how much the non- instantaneous sector throws out to the 
/-channel of the instantaneous sector at each step of the RG evolution. To find out, one 
should make the partial- wave expansion of the function Af in Eq. (|T0|). Since Af(6) is 
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concentrated on values of 9 near 0, the partial-wave coefficients almost do not depend on I, 
provided that / is not parametric on g. Indeed, 



since the Pi are normalized so that Pi(l) = 1. Therefore, the RG equation for f) is identical 
to that for f . The initial value for fi is also independent of I, since it comes from the 
partial-wave expansion of the function ft Tee {9) in Eq. (H) which also peaks near 9 = 0. 
Therefore, to leading order, the RG does not discriminate between channels with different 
angular momenta. If a condensate with nonzero angular momentum is formed, the rotational 
symmetry is broken, like in the A phase of He 3 . 

However, at any finite g, the coincidence of the RG evolution of fi with different I is not 
exact. Moreover, the RG approach does not give us the value of the gap, or the energy of the 
ground state, but merely yields the typical scale of the gap. Therefore, the two gaps with 
different I may have the same asymptotic behavior, but with different numerical coefficients, 
and the corresponding superconducting states may have different energies. At this stage, 
the natural assumption seems to be that the state with I = is favored, and the ground 
state does not break rotational symmetry, but the question of forming a condensate with 
nonzero angular momentum of the Cooper pair should be investigated in a more careful 
manner. We defer this question to future work. We have also left out the possibility of a 
numerical estimation of the gap at moderate densities, which is a very interesting question 
from the phenomenological point of view. Presumably, a reasonable estimation could be 
found by solving the Eliashberg equation of the type described in Appendix [B|. Nor did we 
try to solve the problem at finite temperatures [|]. 
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APPENDIX A: NON-RENORMALIZATION OF THE NON-INSTANTANEOUS 

INTERACTION 

Here we present an explicit check that the soft sector of the theory is not renormal- 
ized during the RG flow. Let us consider one-loop corrections to the scattering amplitude 
p, — p — > p, —p. There are two Feynman diagrams to be considered: 
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The contributions from the two graphs are of the same order, let us evaluate only the first 
graph, which is of order 



g 4 dk dk - TT^—n 2 |fcoh 2 ( A1 ) 

Notice that the momentum of the internal fermion lines should be inside the shell e _1 5 < 
l e p+q| < 8. Both gluon lines are supposed to belong to the non- instantaneous sector, so 
k ^ m^S 1 / 3 . The integral over ko is dominated by k^ < 5. Therefore, the integral in Eq. 



flAl| ) is of order 

Therefore, the correction is of order g A q$ 2 - 

Now assume that the external momenta of the final particles are slightly different from 
those of the initial particles, and the difference is q. The tree-level amplitude is of order 
g 2 q~ 2 ■ During each step of the RG evolution, this amplitude receives a correction of order 
g 4 qs~ 2 , where S is the moving RG scale. During the part of the RG evolution when the 
interaction is non-instantaneous, q$ > q, and one sees that all the accumulated correction 
(~ g 4 q~ 2 ) is still smaller than the tree amplitude by a factor of g 2 . One concludes therefore 
that there is no substantial renormalization of the non-instantaneous interaction during the 
RG evolution. 



APPENDIX B: THE ELIASHBERG THEORY 



The Eliashberg equations |L5] can be considered as the generalization of the BCS gap 
equation to the case of a non-local interaction. Here we will be trying to reproduce only the 
leading exponential behavior exp(— ^=-) of the gap, but not the power (g~ 5 ) part. Presum- 
ably, a more careful treatment of the Eliashberg equation should reproduce the subleading 
g~ 5 factor and give an estimate for the numerical coefficient. 

The generalization of the gap equation of the type written in Ref. JF] to the case of the 
non-local magnetic interaction is0 

a / s 1 A(gg) 

^ P0) 3 9 J (2vr)4 | p _ q |2 + ^JgzgL ■ 4 + e 2 + A 2 (go) ^) 



3 The coefficient § in Eq. (|Bl]) can be understood as follows. If one replaces the gluon propagator 
in Eq. (|Bl|) by i , Eq. ( pi] ) gives the same gap as Eq. (^) if the numerical coefficient in the RHS 
is |. The coefficient one should put in the RHS of Eq. flBl| ) should be twice smaller than that, due 
to the screening of the electric field, which reduces the effective coupling by a factor of 2. For the 
color-flavor locking scheme ]7j, there should be two equation for two gaps; in weak coupling they 
can be written as one equation (|B1| ) which is valid to the leading order. 
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Notice that, as in the usual Eliashberg theory, the gap is a function of po only but not of 
p. Indeed, the dependence of the RHS of Eq. ( [B"T1) on p is only in the gluon propagator, and 
if p is near the Fermi surface any change of p can be compensated by a rotation of q. The 
conventional Eliashberg equations are actually a set of equations for A(p ) and Z(p ), where 
Z(po) is the wave function renormalization, but as explained in Sec. [V], Z(po) « 1. Notice 
also that in Eq. (|B1|) we write /j, 2 instead of | m^, the reason is that we will be working with 
exponential scales so the difference between these two coefficients will not affect the leading 
exponent. Integrating over q, one finds, 

A(p„) = JL J dm ln . m (B2) 



where only the leading logarithm is written. To the leading log, the integral in Eq. can 



be split into two regions, < go < Po and po < qo < /i, where in the first ln | po ^ go | ~ hi ^ 
and in the second ln , — ^ — r « ln Introducing the logarithmic scales 

x = ln— , y = In— , x = ln-^- 
Po Qo A 

where A = A(p = 0), the Eliashberg equation becomes, 

A(x) = ^_(x JdyA(y) + JdyyA(y)^ (B3) 

x 

Differentiating Eq. (|B3|) with respect to x twice, one finds, 



A"(x) = -£- 2 A(x) (B4) 

As the boundary conditions, from Eq. (|B3|) one can check that A(0) = and A'(xo) = 0. 
The solution to Eq. (|B4|) is, 



A(x) = Ansinf — %^x ) 



To satisfy the boundary condition at xo, one requires j^r^o = f , from which one finds 

x ° = v^' anc ^ ^ e S a P a ^ sma ^ energies is A ~ e~ x ° ~ exp(— -^), reproducing the leading 
exponential behavior of our RG result. This is the minimal energy cost to create a fermion 
excitation. The energy- dependent gap A(p ) is 



A(p ) = A sin( — |=-m — 
\3v27r po 



for po > A ~ exp(-^g) 
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